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A new purification scheme is proposed which applies to ar- 
bitrary dimensional bipartite quantum systems. It is based on 
the repeated application of a special class of nonlinear quan- 
tum maps and a single, local unitary operation. This special 
class of nonlinear quantum maps is generated in a natural way 
by a hermitian generalized XOR-gate. The proposed purifica- 
tion scheme offers two major advantages, namely it does not 
require local depolarization operations at each step of the pu- 
rification procedure and it purifies more efficiently than other 
know purification schemes. 

PACS numbers: 03.67.-a,03.67.Hk,03.65.Bz 



I. INTRODUCTION 

Quantum information processing was proven to be su- 
perior to classical information processing in several re- 
spects. The essence of quantum information is its ability 
to employ the linearity of quantum mechanics on compos- 
ite systems for practical purposes. Many of the already 
demonstrated procedures rely on the use of highly 
entangled quantum states. Entangled states are never 
generated in ideal form. Typically, either the source pro- 
ducing entangled quantum states or the communication 
channel with which entanglement is transfered to remote 
parties adds noise. In order to be able to exploit entangle- 
ment efficiently it is desirable to remove as much of this 
additional noise as possible. This goal can be achieved 
by purifying or concentrating entanglement. 

The manipulation of quantum states is realized using 
quantum networks. Usually, they are constructed out of 
simple elements, so called quantum gates. Among them 
the two-particle quantum XOR-gate plays a funda- 
mental role. In this 2-qubit gate, the first qubit controls 
the target qubit: if the control is in state |0), the target is 
left unchanged, but if the control qubit is in state |I) the 
target's basis states are ffipped. Together with one-qubit 
operations this gate forms a universal set of quantum 
gates allowing the implementation of arbitrary unitary 
operations acting on qubits [Q . It has been demonstrated 
that the quantum XOR-gate can be used for many prac- 
tical tasks of quantum information processing, such as 
quantum state swapping Q, entangling quantum states 
P , performing Bell measurements , dense coding [|| , 
and teleportation Furthermore, in combination with 
selective measurements it can be used for implementing 
nonlinear quantum transformations which may be used 
for optimal state identification and for state purification 
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For many practical tasks of quantum information pro- 
cessing it is desirable to extend the basic notion of such 
a quantum XOR-operation to higher dimensional Hilbert 
spaces. Indeed, most of the physical systems that have 
been proposed to hold qubits, such as multilevel atoms 
or ions and multipath-interferometers |l^, could 
equally well encode larger alphabets. However, there is 
a considerable degree of freedom involved in such a gen- 
eralization. 

In the present paper we use a hermitian generaliza- 
tion of the quantum XOR-gate which applies to arbitrary 
dimensional Hilbert spaces and which allows to imple- 
ment a special class of nonlinear quantum transforma- 
tions in a natural way. These nonlinear transformations 
can be used for the preparation of quantum states and 
for efficient quantum state purification. This will be ex- 
emplified by discussing state purification of generalized 
Bell-states. These latter quantum states are of consider- 
able interest in quantum information processing in higher 
dimensional Hilbert spaces. Compared with the other 
known purification scheme which is valid in arbitrary di- 
mensional Hilbert spaces and which has been proposed 
by Horodecki et al. |Q our proposed purification proce- 
dure offers two advantages. Firstly, it does not involve 
a depolarization operation at each step of the iteration 
procedure. Typically, it is not easy to implement such a 
depolarization operation. Furthermore, such a depolar- 
ization operation requires a minimal size of the ensemble 
of quantum states it is applied to. Secondly, it will be 
demonstrated that our newly proposed method is more 
efficient. 

The paper is organized as follows. In Sec. II the hermi- 
tian generalized XOR-gate is introduced which allows to 
implement the special class of nonlinear quantum trans- 
formations needed for our proposed purification scheme. 
The resulting class of nonlinear quantum maps is dis- 
cussed in Sec. HI. The new quantum state purification 
scheme, its basic properties and its efficiency are exem- 
plified in Sec. IV. 



II. A HERMITIAN GENERALIZED GXOR-GATE 



Let us start by summarizing characteristic properties 
of the XOR-gate as they are known for qubit systems. 
For qubits the action of the quantum XOR-gate onto a 
chosen set of basis states {\i)} with i e {0,1} of the 
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Hilbert space of each qubit is defined by 

XORu\^),\J), = \^),\^(BJ), . (1) 

This transformation has the following characteristic 
properties: (i) it is unitary and thus reversible, (ii) it is 
hermitian and (iii) i(Bj = if and only ifi = j. The first 
(second) index denotes the state of the control (target) 
qubit and ® denotes addition modulo(2). 

Let us now consider the problem of generalizing the 
quantum XOR-gate to higher dimensional Hilbert spaces. 
The desired generalized quantum XOR-gate (GXOR- 
gate) should act on two D -dimensional quantum sys- 
tems. In analogy with qubits we will call these two sys- 
tems qudits. The basis states \i) of each qudit are labeled 
by elements in the ring Zjj which we denote by the num- 
bers i — 0, I? — 1 with the usual rules for addition and 
multiplication modulo(D). In principle, the GXOR-gate 
could be defined in a straightforward way by using Eq. (|^) 
and by performing i® j modulo(D), i. e. 

GXORu\i),\j)2 = \^)i\^®j)2- (2) 

The GXOR-gate defined in (|^) is unitary but not her- 
mitian for D > 2. Therefore it is no longer its own 
inverse. Thus, the inverse GXOR-gate has to be ob- 
tained from the GXOR-gate of Eq.(||) by iteration, i.e. 
GXOR^2^ = (GX0i?i2)^~i = GXOrI^ ^ GXOR12. 
All these inconvenient properties of this preliminary def- 
inition (H) can be removed by the alternative definition 

GXOR,2V),\])2 = \^iVQ3)2- (3) 

In Eq.(|[) iQj denotes the difi'erence i — j modulo (D). In 
the special case of qubits the definition of Eq. (^) reduces 
to Eq.(|l|) asiQj = i(Bj modulo{2). Furthermore, this 
definition preserves all the properties of Eq.(0) also for 
arbitrary values of namely it is unitary, hermitian and 
i Q j = modulo{D) if and only if i = j. 

The GXOR-gate of Eq. admits a natural extension 
to control and target systems with continuous spectra. 
In this case the basis states \i) are replaced by the basis 
states {|a;)} with the continuous variable x € Tl. These 
new basis states are assumed to satisfy the orthogonality 
condition {x\y) = 5 (x — y). Furthermore, as the dimen- 
sion D tends to infinity the modulo operation entering 
Eq.(||) can be omitted. Thus, for continuous variables 
the action of the GXOR-gate becomes 

GX0RMi\y)2^\^)i\^-y)2- (4) 

Let us note that this definition for the case of continu- 
ous variables is different from the generalized XOR-gate 
proposed in Rcf. This latter gate is not hermitian 

whereas the GXOR-gate of Eq.(Q) is both unitary and 
hermitian. The GXOR-gate of Eq. (g) can be repre- 
sented in terms of a translation and a space inversion, 
namely 

GX0R,2 \x), = n2e(-^^^''**'') \x), \y)^ . (5) 



Thereby Py denotes the canonical momentum operator 
which is conjugate to the position operator y^^^ acting on 
quantum system 2 and II2 is the corresponding operator 
of space inversion. 

Let us discuss a possible physical realization of the 
GXOR-gate defined by Eq.(^) which is based on nonlin- 
ear optical elements. For this purpose we assume that 
the two quantum systems which are going to be entan- 
gled are two modes of the radiation field. The basis 
states {i = 0,...,D — 1) of the first quantum sys- 
tem are formed by n-photon states of mode one with 
< n < D — 1. The basis states of the second quantum 
system |fc)2(fc = 0,...,D — 1) are formed by Fourier trans- 
formed n-photon states of this latter mode, i.e. \k)2 = 

exp{i2Trkn/D)\n)2- Let us further assume 
that the dynamics of these two modes of the electromag- 
netic field are governed by the Kerr-effect Thus, 
in the interaction picture their Hamiltonian is given by 
H — hx<i\cLia\a2 with the creation and annihilation op- 
erators a\ 2 and ai_2 of modes 1 and 2, respectively. For 
the sake of simplicity the nonlinear susceptibility x is as- 
sumed to be real- valued and positive. Preparing intitially 
both quantum systems in state |i)i|fc)2 after an interac- 
tion time of magnitude t — 2t:/{Dx) this two- mode sys- 
tem ends up in state |V')i2 = \i)i\k — i)2- Applying to this 
latter state a time reversal transformation which may be 
implemented by the process of phase conjugation n& we 
finally arrive at the desired state — k)2- Thus this 
combination of a Kerr-interaction with a time reversal 
transformation is capable of realizing the GXOR-gate of 
Eq. d). 

As the GXOR-gate of Eq.(|) differs from the alter- 
native definition of Eq.(^) by local unitary operations 
only both are expected to exhibit similar properties as 
far as entanglement operations are concerned. However, 
with the help of this GXOR-gate alone already a variety 
of interesting quantum operations can be implemented 
without having to use additional local unitary transfor- 
mations. In addition, as will be shown in Sees. Ill and IV 
this entanglement operation is particularly useful in the 
context of quantum state purification for implementing 
nonlinear quantum state transformations. 

As a first application let us consider the preparation 
of a basis of entangled states from separable ones. If 
|^)|m) with l,m,= 0,...,D — 1 denotes an orthonormal 
basis of factorized states an associated basis of entangled 
two-particle states is given by 

\i:i^) = GXORi2mi))i\m)2]. (6) 

Thereby F denotes the discrete Fourier transformation, 

D-l 

i.e. F\l) = [I/VD) J2 exp{i2TTlk / D)\k) . For qubits this 

fc=0 

unitary quantum transformation leads to the well known 
basis of four Bell states. In the simplest higher dimen- 
sional case of I? = 3, for example, the first four states of 
this entangled generalized Bell basis are given by 
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-^[|00) + |11) + |22)], 
v3 

^[|00) +e-^2'^/=^|ll> +e'27r/3|22)]^ 
v3 

-i=[|02) + |10) + |21)]. 



(7) 



As the GXOR-gate is hermitian it can also be used to dis- 
entangle this basis of generalized Bell states again by in- 
verting Eq.(^. This basic disentanglement property is of 
practical significance. It enables to reduce Bell measure- 
ments to measurements of factorized states. Examples 
where these latter types of measurements are of central 
interest are dense coding |^] and quantum teleportation 
schemes 

The basis of entangled Bell states resulting from Eq. (|^) 
can be used for teleporting an arbitrary D-dimensional 
quantum state from A (Alice) to B (Bob). For this pur- 
pose let us assume that A and B share an entangled 
pair of particles prepared in state |V'/m) as defined by 
Eq. (^ . If A wants to teleport an unknown quantum state 

D-l 

\x) = Q;„|n) to B she has to perform a Bell measure- 

mcnt which yields one of the entangled basis states of 
Eq.(^ as an output state (compare with Fig. (|l|)). Con- 
ditioned on the measurement result of Alice, Bob has to 
perform an appropriate unitary transformation onto his 
particle which prepares this latter particle in state \x)- 
This arbitrary dimensional teleportation scheme rests on 
the identity 



D-l 



-i27Tjm/ D 
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Ui^\n) = e-''2™('-j)/-D|„ - fc - m). 



(8) 



This basic relation for teleportation for an arbitrary di- 
mensional state Iy) can be derived in a straightforward 
way from Eqs. (^ and (^. The classical communication 
requires 21og2(-D) bits, which is the minimum necessary 
in all quantum teleportation schemes. 



III. NONLINEAR QUANTUM MAPS ON 
DENSITY MATRICES 

With the help of the hermitian GXOR-gate of Eq.(|) 
an interesting class of nonlinear quantum maps can be 
implemented in a natural way. Together with filtering 
measurements acting on a target quantum system t the 
GXOR-gate of Eq. (|3[) induces nonlinear transformations 
of quantum states of a control system c. This can be 
demonstrated most easily by considering the case of two 
qudits which are prepared in the quantum states tr* and 
a'^ initially. Let us perform the quantum operation 



^(c^^^7* 



Tr[A (cr'= (g) cr*) At] 
on these two qudits with 

A = (Ic ® P) GXORct. 



(9) 



(10) 



Thereby Ic denotes the identity operator acting in the 
Hilbert space of the control system and P = \p)^^ {p\ is 
the projector onto the state \p)^ of the target qudit. With 
the decomposition 
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Eqs. (p) and (10) may be rewritten in the form 
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E '^tAi \^)cc{j\®P\iek),,{jei\P 



ikl 



(12) 

Assuming that both control and target qudit are pre- 
pared in the same state initially, i.e. a'^ = cr*, it turns 
out that Eq.(^) is equivalent to the relations 

^(cT^a*^a^)-<„,p„,®P, 

D-l 

E <j<-p.i-pN>cc(jl 

• (13) 



^ output D — 1 

i 

As a result of the quantum operation (|^) the combined 
system formed by the control and the target qudit forms 
a factorizable state with the target qudit being in state 
|p)(p|. According to Eq.dl^) the density matrix ele- 
ments of fj'^ with respect to the computational basis |i) 
[i = 0, D ~ 1) have been multiplied with each other. 
The final state of the control qudit is prepared with prob- 

D-l 

ability = E '^i i^i~p i~p- 1^ l^ie case of qubits, for 

i 

instance, the nonlinear transform of the (unnormalized) 
control density matrix is given by 



output 



(<^li)' 



(^lo)' 



(^22)' 



(14) 



if projected onto the state |0)t. The elements of the orig- 
inal density matrix have been squared. If one projects 
onto the state |l)t, the original density matrix elements 
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are mixed in a more complicated way and one obtains 
the (unnormalized) density matrix 



output 



33 



fJinCr: 
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(TonCr 



V 



20^^13 



'^01 "'32 

<72lO-12 



^02 "'31 
f^l2 0-21 
^12 0-11 



30"03 "31'^02 '^32^01 "33"00 



'''03 '''30 
1^13 '^20 
1^23 '''10 



(15) 



From Eq. ([1^) it is easy to verify that the quantum 
operation (H) has the following basic properties: (i) it 
maps density matrices onto density matrices, (ii) it is not 
injective and nonlinear, and (iii) there are states invariant 
under this transformation. It is also possible to extended 
the quantum operation of Eq. (^ to cases in which there 
is more than one control system and in which both the 
control and the target systems are composite quantum 
systems each of which consists of M qudits. In this case 
a'^ describes a general M-qudit state of the form 



(16) 



with i = and j = (ji, ja/)- In Eq.(||) the 

operator A has to be replaced by 



A = (Ic P)nfUllf^^GXORl: 



(17) 



with the projection operators P — YiiLi ®Pti and Pt- = 
|pi)j.(. (pi| onto state |Pi)j. of the M-qudit target system 

U. Thereby the GXOR-gate GXOR^^l operates on the 
j-th qudit of the control and of the i-th target system. 
For the special case of = for all ti, for example, the 
resulting final state of the control system is given by 



EH) 

ij 



c \l+N I 



output 



(18) 



i 

and is prepared with probability Pc = X^(o'ii)^^^- 



IV. BIPARTITE PURIFICATION IN HIGHER 
DIMENSIONAL SPACES 

In general, for A'^ = 1 the nonlinear quantum trans- 
formation of Eq. (|l8|) has not only invariant states for 
Pti = |0)titi(0| but also for other projectors Pt^. This 
suggests to use this nonlinear quantum transformation 
for the purification of quantum states of a two-qudit sys- 
tem. For the special case of a control system consisting 
of two-qubits such a purification scheme which is based 
on the nonlinear quantum transformation of Eq. ( p^ has 
already been proposed previously [ll| . In order to dis- 
cuss an analogous purification scheme in arbitrary dimen- 
sional Hilbert spaces we start from the observation that 
for M = 2 the entangled basis state IV'oo) of Eq.(||) is a 
fixed point of the nonlinear two-particle quantum map 



of En . (| i 8| ) . Inspection of the two 'squared' matrices (|TJ) 
and ( 15]) shows that this is true not only for the projector 
Pt^ = |0)tt(0| but also for Pt, = |l)tt(l|. Therefore, this 
nonlinear quantum transformation may be used for puri- 
fying quantum states towards the entangled state I'i/'oo)- 
Thereby the possibility to use both projectors in the pu- 
rification process increases its efficiency considerably. 

Purification in higher dimensional Hilbert spaces has 
been considered previously by Horodecki et al. Q. 
These authors generalized the approach of Bennett et 
al. to arbitrary dimensional Hilbert spaces. The 
purpose of their protocol is to distill Bell states from a 
noisy channel. Their protocol combines two basic steps, 
namely a nonlinear quantum map which 'squares' the 
density matrix elements and a depolarizing channel con- 
verting the resulting output state into a Werner state. 
The depolarizing channel guarantees that at each step 
of the purification protocol an initially prepared Werner 
state is mapped again onto a Werner state with a higher 
admixture of the Bell state. This protocol is capable of 
purifing all non-separable Werner states in arbitrary di- 
mensional Hilbert spaces. The depolarization operation 
involved in this purification protocol requires that suit- 
ably chosen local unitary operations have to be applied 
to a sufficiently large number of two-qudit systems. For 
the case of qubits it has been shown that such a depolar- 
ization may be achieved with a set of 12 suitably chosen 
unitary operations ||l^. For qudits with D > 2 suitable 
minimal numbers of unitary operations are not known at 
present. In view of these inconveniences in implementing 
a depolarizing operation it appears desirable to develop 
alternative purification strategies which do not involve 
such a depolarization procedure. For the case of qubits 
such a procedure has already been developed by Deutsch 
et al. In the following we propose such a method 

which applies to arbitrary dimensional Hilbert spaces. 
It is based on the nonlinear quantum transformation of 
Eq.(|l8|) (with = 1) followed by a single local unitary 
operation acting on both qudits. Thus, this proposed 
purification procedure does not require a minimal num- 
ber of two-qudit systems at each step of the purification 
protocol. 

In order to exemplify basic properties of our purifica- 
tion scheme let us consider the purification of a Werner 
state of the form 



A|Voo)(^oo| + (l-A)l/2?' 



(19) 



where the parameter A is related to the fidelity F 
("000 1 c*^ I'i/'oo) through the expression 



F = X 



1- A 
£|2 ■ 



(20) 



The state (|l|) may result from a physical situation where 
two spatially separated parties, say A(lice) and B(ob), 
want to share the entangled basis state I'i/'oo) but with 
a probability of (1 — A) the transmission of this entan- 
gled pair through a quantum channel leads to unwanted 
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noise represented by the chaotic state This ini- 

tial quantum state a'^ is non-separable if and only if 
A > \d — (1 + -D)^^ so that a purification scheme 
can succeed only for those values of A. 

We propose a purification scheme which is based on 
a two-step iteration procedure. In the first part of each 
iteration step we apply the nonlinear quantum map of 
Eq.(|l8|) (with iV = 1) by projecting onto an arbitrary 
two-qudit target state (i = 0,...,D — 1). Corre- 

spondingly, the initially prepared Werner state ( |l9| ) is 
converted into the quantum state 



output 



Ai 1-000) (^00 1 



D-l 

A2l/i52^A3^ |^^)(n| (21) 

1=0 



with the coefficients A; depending on the initial choice 
of A. However, it turns out that both the pure state 

D-l 

IV'oo) (V'oo I and the mixed state ^ are fixed points 

1=0 



of the nonlinear quantum transformation of Eq.(18). 
Therefore, an additional local unitary transformation is 
required in order to guarantee convergence of the itera- 
tion procedure towards the desired final state |'(/'oo)(V'oo|- 
Thus, in the second part of each iteration step parties A 
and B perform a local twirling operation ||l^ Ua ®Ug, 



I.e. 



c(2) 



Ua' 



U*bCjI%^,U\®U*1 (22) 



Ultimately this transformation achieves convergence to- 
wards our desired final state |'0oo)('!/'oo| by altering 
the mixed state but still leaving state 

|'(/'oo)('0oo| invariant. Therefore, the depolarization op- 
eration of the protocol of Horodecki et al. |Q is replaced 
by a single twirling operation. For the unitary transfor- 
mation involved in this twirling operation we propose to 
choose a discrete Fourier transform, i.e. 



UA{B)\k)A{l 



1 



D-l 



n=0 



exp 



(i2Tvkn/D) 



n)A(B) 



(23) 



This choice is motivated by the desire to increase the suc- 
cess probability of the purification process and to maxi- 
mize the radius of convergence of the iteration procedure. 
As both the intermediate output state of Eq.(21) and 
the unitary transformation of Eq.(p3|) are invariant under 
transformations of the basis states of the form \i) |»+1) 
{i = 0, ...,D — 1), all the projections Pt^ = yield 
the same success probability thus increasing efficiency. 
Iterating the two-step procedure based on the nonlinear 
quantum transformation of Eq. ( p^ and the local twirling 
operation of Eq. (^) yields our proposed purification pro- 
cedure. 

Numerical simulations performed for dimensions 2 < 
D < 20 demonstrate that this purification procedure 
is capable of purifying almost all non-separable Werner 
states of the form of Eq. ( 19 ) . The dependence of its range 



of convergence on the dimension D of the Hilbert space 
is apparent from Fig. 3. The dashed line indicates the 
minimum values Fc — 1/-D of the fidelity for which the 
Werner state of Eq.(^9|) is still non-separable. The pu- 
rification protocol of Horodecki ct al. jl^ converges for 
all initial values F > Fc- The solid line represents the 
minimal initial value of the fidelity for which our protocol 
purifies. It is clear from Fig. 3 that the range of conver- 
gence of our purification scheme is slightly smaller but 
approaches the ideal limit Fc with increasing dimension 
D of the Hilbert space. Already for moderately large 
dimensions D our range of convergence approaches the 
ideal range closely. 

Let us now compare the efficiency of our purification 
protocol with the one proposed by Horodecki et al. [0 . 
In particular, we are interested in answering the question, 
how many iterations are needed to obtain state {ipoo) with 
a prescribed final fidehty F final = (V'ook/„a; I'^oo) for a 
given dimension D7 In order to clarify the calculation 
of this efficiency let us briefly reconsider the basic steps 
involved in a purification protocol. They are represented 
schematically in Fig. 4. Initially the purification pro- 
cess starts with an ensemble of Niniuai Werner states 
each of which is described by Eq.(p9[). In each step of 
the iteration procedure the ensemble of two-qudit states 
is divided into two equal parts which serve as control 
and target systems. The nonlinear quantum transforma- 
tion is performed by projecting onto one of the target 
states = \n){ii\ [i = 0,...,D — 1). As our initial 
state and our unitary transformation of Eq.(^3|) arc in- 
variant under the transformation \i) — > |i-t-l) of the basis 
states, all these projections are equally probable despite 
the fact that our purification procedure does not yield a 
Werner state at each step of the iteration procedure. This 
nonlinear quantum transformation is followed by a local 
twirling transformation (compare with Eq.(p^)) which is 
based on the discrete Fourier transformation of Eq. (^3|) . 
Thus, after the first step of our purification procedure we 
are left with [pi x Ninitiai/'2\ purified two-qudit systems. 
Thereby pi denotes the probability of obtaining the tar- 
get qudit in one of its basis states \i) {i = 0, D — 1). 
Continuing this iteration procedure after n iterations the 
number of remaining purified two-qudit systems is given 

by 



N, 



purified 



initial 
2n 



1=1 



(24) 



Accordingly, the efficiency 77 of this purification process 
is given by 77 = Npurified/Ninuiai- 

In Fig. 5 and Fig. 6 the dependence of the efficiency 
ij on the initial fidelity of the Werner state F (com- 
pare with Eq.(pO|)) is depicted for dimensions D = 6 
and D = 9 and for different values of the final fidelity 
Fftnai = ('0oo|o-^i„a/lV'oo)- From these figures it is ap- 
parent that our protocol requires fewer steps than the 
protocol of Horodecki et al. |jl^. Furthermore, numerical 
studies also indicate that in both purification protocols 
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the success probabilities pi entering Eq.(p^ are compa- 
rable in magnitude. Thus, the overall better efficiency of 
our purification protocol which is apparent from Figs. 5 
and 6 reflects the fewer number of steps n which are re- 
quired for achieving a given final accuracy. With increas- 
ing accuracy of the final purified state this difference in 
efficiencies between both purification protocols becomes 
larger and larger. From Figs. 5 and 6 one also notices a 
second characteristic feature which has been found also in 
other numerical simulations. For a given value of the fi- 
nal fidelity the differences between the efficiencies of both 
protocols becomes smaller with increasing dimension D 
of the Hilbert space involved. 

V. CONCLUSIONS 

A novel purification scheme has been proposed. It is 
based on the iterative application of a special class of non- 
linear quantum maps and a single, local unitary trans- 
formation. The required nonlinear quantum map can 
be implemented conveniently by a hermitian generalized 
quantum XOR-gate. The proposed purification scheme 
has several attractive features. Firstly, it applies to arbi- 
trary dimensional bipartite quantum systems. Secondly, 
it does not require a depolarization operation at each step 
of the iteration procedure. It rests on a single twirling 
operation which is performed at each step of the iterative 
purification scheme. Thirdly, the proposed procedure 
achieves purification in a very efficient way. In particu- 
lar, it has been demonstrated that it achieves purification 
of Werner states in a more efficient way than the other 
known purification protocol which has been introduced 
by Horodecki et al. |jlj]. Furthermore, its almost maxi- 
mal range of convergence indicates that the employed lo- 
cal twirling operation which is based on a discrete Fourier 
transform is a good choice. Such a transformation can 
be implemented easily in many quantum systems. 

The proposed purification method may also be gener- 
alized to multi-partite quantum systems. In this context 
it would be particularly interesting to develop efficient 
purification protocols for GHZ-like quantum states. 
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FIG. 1. Schematic representation of 
scheme involving Bell measurements onto the generalized Bell 
states of Eq.(^. 
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FIG. 2. Schematic representation of the GXOR-gates and 
projections involved in the nonlinear quantum transformation 
of Eq. dlil). The qudits are represented by dots. The dots 
of the first line represent the M qudits of the control system. 
The dots of the following lines represent the M x A*' qudits of 
the target systems ti,t2, In- The GXOR-gate GXOB!~Jt\ 
acts on the j — th qudit of the control and target system ti 
with j G {1, 2, A'/} and i € {1, 2, N}. The operator P- 
projects the state of the whole systems onto state |0) (0| with 
|0> = |0>JO)2...|0)„^. 
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FIG. 4. This scheme illustrates the calculation of the ef- 
ficiency r]. The initial collection of qudits is split into two 
equal parts. One part forms the control, the other the target 
qudits. The success probability of the process after the first 
step is pi, the total number of 'squared' qudits is given by 
{l/2)Ni„itiaiPi- The procedure is repeated until the required 
fidelity Ffi„ai for the resulting state is reached after n steps. 



FIG. 3. Dependence of the minimal initial fidelity Fc 
needed to purify a Werner state (compare with Eq. (p^) 
as a function of the dimension D: our protocol (full line), 
Horodecki's protocol (dashed line). 




FIG. 5. Dependence of the efficiency -q (compare with 
Eq. (p^) on the initial fidelity F for a fixed final fidelity 
Ffinai = 1 — 10"'' and for dimension D = 6. The solid line 
gives the results of the proposed method and the dashed line 
the results of Horodecki's protocol. 
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FIG. 6. Dependence of the efficiency rj (compare witfi 
Eq. (p^)) on the initial fidelity F for a fixed final fidelity 
Ffinai = 1 — 10~^ and for dimension D = 9. The solid line 
gives the results of the proposed method and the dashed line 
the results of Horodecki's protocol. 
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